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O ! Abstract 

o. 

We have shown that the B — L generation due to the decay of the therniahy 
produced superheavy fields can explain the Baryon assymmetry in the universe if 
the superheavy fields are heavier than 10^'^"^^ GeV. Note that although the super- 
' r~| , heavy fields have non- vanishing charges under the standard model gauge interactions, 

C^l the thermally prduced baryon asymmetry is sizable. The B — L violating effective 

pi i' operators induced by integrating the superheavy fields have dimension 7, while the 

operator in the famous leptogenesis has dimension 5. Therefore, the constraints from 
the nucleon stability can be easily satisfied. 
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1 Introduction 

To understand the origin of Baryon number B in the universe is one of the most interest- 
ing subjects in the particle cosmology. The abundance of the Baryon in the universe is 
estimated by the nucleosynthesis analysispQ and is observed by the WMAPj^, and it is 
quite impressive that they have given the consistent value for the Baryon density in the 
universe, which is roughly 

^ ~ 10-1°, (1.1) 

s 

where n^ and s are the Baryon number density and the entropy density, respectively. After 

Sakharov[3] pointed out the three conditions for the generation of the Baryon number B in 

the universe, many mechanisms for baryogenesis have been studied in the literature [H [5|, 

EllTj. One of the most attractive scenario for the baryogenesis is the GUT baryogenesis [1] 

in which the decay of superheavy gauge bosons and Higgs appeared in the GUT produces 

the baryon number. Unfortunately, the produced baryon number is known to be washed 

out by the sphaleron process[8] in the standard model (SM). Since the sphaleron process 

conserves the B — L number, it is important to produce non- vanishing B — L number. The 

most famous scenario to produce B — L number is the leptogenesis|5j, in which the lepton 

number L is produced by the decay of the right-handed neutrino. Especially, thermal 

leptogenesis, in which the lepton number is produced by the decay of the right-handed 

neutrino produced thermally, is one of the most interesting scenario because the observed 

baryon number can be related with the measurements on neutrino masses and mixings. 

However, since the scenarios in which the thermal leptogenesis can be applied are limited, 

other possibilities to produce the B — L number are worth considering. In this paper, we 

study the B — L production by the decay of certain superheavy fields with intermediate 

masses, which can be a remnant of some GUT models. 



2 B — L Violating Interactions 

In the SM, the renormalizable operators cannot break the B and L numbers. Therefore, 
the non-conserving interactions appear in the higher dimensional operators whose mass 
dimensionis larger than four. For example, the dimension five operators llhohn between 
the doublet lepton / and the doublet Higgs /i^ have non- vanishing B — L charges and give 
neutrinos masses. The dimension 6 operators qqu^^e^^ break the B and L numbers, which 
can induce the proton decay. (For our notation of the particle contents, see Table [H) 
llhohr) violates L and B — L, while qqu'^^Cn violates B and L but not B — L. These higher 
dimensional operators, Uhnho and qqu'^e^, can be induced by integrating the superheavy 
right-handed neutrino v'^ and the GUT gauge boson X, respectively, as in FigUl The 
right-handed neutrino plays an important role in the leptogenesis scenario. And the X 
gauge boson also plays a crucial role in the GUT baryogenesis. Therefore, in order to 
produce B — L number, it must be important to understand which superheavy particles 
can induce the B — L violating higher dimensional operators. 

Now, we discuss on other B — L violating operators than llh£,ho- In the literature, in 
the context of the nucleon decay, B and/or L violating operators have been classified in 
the SM^ and in the minimal supersymmetric SM (MSSM)p^. In the SM, there is no 
dimension six i? — L non- conserving operator. It is in dimension seven that we can find 
out B — L non-conserving operators 
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Table 1: The particle contents and charges. 
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Figure 1: The decomposition of UhDhc (upside) and qqu^j^e^^ (downside). The fook over 
operators means the contraction. 



qdytlho , u'}id'}id'y,hD , e'^^lllhD , qqd'^^t^ Wj~, , quytW^h}j~, , qe'^^d'^^d'^^hy 



D ; 



'^r^r^r'- "-d 



dj^dj^dj^lhjy , qd^d^l 



u^d^l I 



e" d^'^d'^^d''^ 
^R^R"'R"'R ■ 



The last three operators include a differential operator or a gauge field. Since the differen- 
tial operator can be replaced by the light fermion mass by using the equation of motion, 
the contribution of these operators become negligible and we do not consider the last three 
operators in the foUowings. 

Which particles can induce these B — L violating higher dimensional operators? To 
answer this question, let us decompose these operators into two parts. It is useful to write 
down these operators with the SU{J)) complete multiplets, 10 = (g,M^,e^), 5 = {d'j^,l), 
and 5s = {Ht, he), where Ht is a colored Higgs, as 10 ■ 5 • 5 • 5 • 5^, 10 • 10 • 5''' • 5''' • 
5 J, and 5 • 5 ■ 5 ■ 5 ■ 5j. First of all, supposing that the superheavy fields are scalar. 
Then each part must includes two fermions. Therefore, the decomposition is limited. For 
example, the operator 10 ■ 5 ■ 5 ■ 5 • 5^ can be docomposed as [10 ■ 5 ■ 5^ + 5 ■ 5] or [10 ■ 
5 + 5 ■ 5 • 5s]. For simplicity, we assume that the superheavy fields are included in the 
SU{5) multiplets, 1, 24, 10, 5. (Though it is straightforward to extend the superheavy 
fields with the general representations, we do not discuss the extension in detail in this 
paper.) For the former decomposition, the superheavy scalar belongs to 10 representation 
of SU{5), and for the latter, 5. The concrete decompositions of the operators, qdj^^llhu 
and u'}^d'}^d'y,hD, can be seen in Figl2] and [3l respectively. We denote the superheavy 
fields as the large characters of the SM fields which have the same quantum numbers 
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Figure 2: The decompositons of qd'^^llho- 



under the SM gauge interactions. (In FigJS], the superheavy field A has charges of (3, 2)7/5 
under SU{'i)c x SU{2)l x t/(l)y, which belongs to 45 of SU{b).) It is obvious that the 
superheavy scalars whose representations are 10 or 5 also induce the other B — L violating 
operators. Therefore, we can consider the B — L generation by the decay of the superheavy 
scalar fields which belong to 10 and/or 5 of SU{b). We will return to this scenario in the 
next section. 

If the superheavy fields are fermions, these operators must be decomposed as three 
fermions and one fermion. For example, the operator 10 ■ 5 ■ 5 ■ 5 ■ 5^ can be decomposed as 
[10 ■ 5 ■ 5 + 5 ■ 5s] or [5 • 5 • 5 + 10 ■ 5^] . For the former decomposition, the superheavy fermions 
belong to 1 or 24, and for the latter, they belongs to 10 and the complex conjugate. 
Actually the right-handed neutrinos which belongs to 1 of SU{5) can induce some of these 
operators. Though it must be possible to induce non- vanishing B — L number by the decay 
of these superheavy fermions, we do not discuss this possibility in more detail. We may 
return to this subject in futurq^l. 

In this section, we have decomposed the dimension seven B — L violating operators into 
two parts. By this decomposition, we can address the generation of the B — L number by 
the decay of the intermediate superheavy fields. Though some of the operators obtained 
by the decomposition have still higher dimension than four, we do not decompose them 
further because we do not need the origin of the operators to discuss B — L generation. 



3 B — L Number Generation in the Early Universe 

In this section, we study the B — L generation by the decay of the superheavy scalar fields. 

First, let us fix the particle contents. As discussed in the previous section, some ad- 
ditional fields are needed, and we introduce bosons denoted as Q, U^, D"^, E"^, L whose 



^ If the superheavy fields are vector, these operators must be decomposed as two parts which include 
fermion and anti-fermion. Such decomposition is possible for the operator 10 • 10 • 5^ • 5''' • 5j. Actually, 
the vector bosons which belong to 10 of SU{5) can induce the operator. 
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Figure 3: The decompositions oi uj^d'j^d'^jilhD- 



Interaction 
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f/^ 


E"" 


D^ 
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dim. 4 
dim. 5 


+4/3 
-2/3 


+2/3 
-4/3 


+2 



+2/3 
-4/3 



-2 



Table 2: The generated B — L number by decay of additional particles. 



charges are the same as the Standard Model fermions, g, m^, d^, e^, /, respectively. 

Next, we write down all the dimension four and five interactions which include only 
one superheavy scalar as 

• dim. 4 : 

gc?^L, e^/L, g^wj^ , (/i.e.), 

• dim. 5 : 

d^phoU" , qe^h^j^U" , wg/f/il^f/' , qWJih^j^E^ , q^d^h^j^E" 
eph^jyE", dph^jjD'', in^h^^L, {h.c), 

where we omited the indices of spinor, gauge of 5't/(2) and S't/(3) for simplicity. It is 
interesting that the B — L number of the final states by the decay of Q, V^, E"^, D'^, L is 
a fixed value for each superheavy fields and each dimension of the interactions as in Table 
EJ Since the B — L number of the final states induced by the dimension four interactions 
is different from that by dimension five interactions, the decay can produce non-vanishing 
B — L number. 

For the estimation of the B — L number, it is useful to calculate the mean net B — L 
number e; 

ei = J2 ^^~^f [K^ ^ /) - rii ^ /)] , (3.1) 

/ 



J + 



Figure 4: The Feynman diagrams for calculating ej. 



where i is the initial decay particle with mass of rrii, which contains Q, U^, E'^, D^, L. f 
means the decay modes from the decay of i. Xi_^f is B — L number within the decay modes 
/. r is the branching ratio; i or f means CP transformated state, i.e., anti-particles, e^ 
means generated B — L number for the decay of two particles i and i. Therefore, we can 
obtain the B — L number density ub-l from the number density of the i particle rii and 
the ej parameter as ns-L ~ Ci'^i- After the sphaleron process, the B number density ub is 
obtained as 

ub ~ 0.35?7.B_2, ~ 0.35ejnj. (3.2) 

Therefore, in order to obtain the B number in a comoving frame B = ^^^,we have to know 
ej and n,. 

For the calculation of e^, we denote couplings as follows; 



l- 



a (SM) a (SM) 

b (SM) ^^ hn or /it, 

where y and A are dimensionless couplings, and A is the scale of the higher dimensional 
interactions. By calculating the Feynman diagrams in Fig. SI we obtain e^ as 



2 m? m. 



- J2 I"^ (^iafe^L?/]cd2/irfa) ■ / i^V^l) ■ (3-3) 



(87r)3 A2 levrF, 

(See Appendix for the detail calculation. As an example if we take i = Q and j = If^, the 
summation becomes ^uid>^QddyuddyQdi + ^uiu>^QduyuddyQdi- We can also take j = E"", D", L 
for i = Q. Of course we can take f/^, E'^, D'^, L as i. Once we fix the concrete fields as i 
and j, the factor due to the number of freedom in the loop is appearing. We just ignore it 
in the following for simplicity.) Fj is the total decay width; 

rrii 'sr-^ / I \ rrii mf sr-^ ( \\ \ \ 

^^ = I^ Z^ \yiahyiab) + ^(^J2 1^ \\ab^iab) , (3.4) 

a,b a,b 

where the first term is the contribution from the two body decay, and the second term is 
from the three body decay. Function / in (13. 3 p is the loop function as follows; 



/(a) = l + 2a[l-(l + a)ln(l + l/a)] (3.5) 

1 + Oia) (a A) 

3^ + Oil/a') (a^l) ■ ^- > 



If rrii ~ rrij but rrii < rrij, then the function / is rougly (9(0.1). 

Supposing that the only one couphng dominates the others for each i particle and for 
each dimensional operator. Namely, there are four couplings, t/i, yj, Aj, and Xj. Then, the 
eqs. (13.31) and (13. 4p can be rewritten as 



mf rrij 



(87r)3 A2 I67rr, 

rrij , ,o rrij mf 



■Im [x,Xly}y?j ■ f {mym^) , (3.7) 



lovr 3(87r)'^ A^ 

Moreover, if we take y = \yi\ ~ \yj\ and A = |Aj| ~ \Xj\ and the branching ratio of two 
body decay is comparable to that of three body decay, i.e., y ~ miX/{A\/QTcA), then we 
can obtain simpler equations as 

r- ~ 2^^yl (3.9) 

^ y'f sin 5, (3.10) 



IGtt 



where sinS = Im (AjA||/jyj j /(y^A^). 



Next, let us estimate the abundance of the i particle, rii, and the Baryon number in a 
comoving frame B in the following two cases. In the first case, the particle i is thermally 
produced, the freeze out occurs when the particle i is still relativistic, and no entropy is 
produced by the decay (case A). (Therefore, we assume that the reheating temperature due 
to the inflation is larger than the mass of the superheavy particle i. We discuss whether 
the particle i is still relativistic or not at the freeze out in section 4.) Then, Yj = — is 
given by Yj ~ 0.278^^^^, where the entropy density s is obtained as s = "^g^sT^ and Qeff 
is the number of freedom of i. Therefore, we can obtain 

5^0.1^2ZZ!l^3.5xlO-V, (3.11) 

g*s 

where in the last similarity we use eq. (I3.10p . /sin (5 ~ 0.1, g^s ~ 100, and g^jf ~ 6 for 
i = Q. Therefore, 

?/ ~ 2 X 10"^ (3.12) 

is required to obtain B ~ 10^^°. An additional condition Fj > {avi)ni is required so that 
the estimation hb-l ~ ^i^i is valid, where a and Vi are the cross section of the annihilation 
process and the velocity of the particle i, respectively. If {avi) ~ O.OlTo/mf , this condition 
is roughly rewritten as 



Mpi 



> 10~V- (3.13) 



Here, the decay temperature T^ is deflned by the temperature of the universe when the 
age of the universe is around the lifetime of the particle i, which is given by To ~ 
^yTiMpi/l.66y/g:, where Mpi is the Planck mass as Mp^ = 1.22 x 10^^ GeV. From the 
eq. (I3.12J) . the inequality (I3.13P is rewritten as 

mi > 2 X 10-^Mp; ~ 2 X 10^^ GeV. (3.14) 

When rrii ~ 10^^(10^^) GeV, F^ ~ IQ-^mi ~ 10^(10^) GeV, and therefore, T^ ~ 3 x 

6 



10^^(3 X 10^^) GeV. The higher dimensional couphng is given as A/A ~ AyQixy/nii ~ 5 x 
10-16(5 X 10-i8)(GeV)-i. 

As the second case (case B), we consider the situation in which the density of i and i 
fields dominates the density of the universe. Generically, thermal abundance of the heavy 
particle with long lifetime becomes large and sometimes dominates the energy density of 
the universe. Then, 

P = Pi + Pi = ^miUi = ^9*Tr, (3.15) 

where pi, pj, Ui, (?*, and Tr are the energy density of i field, that of i field, the number 
density of i field, the total number of effectively massless degrees of freedom, and the 
temperature after the i and i field decay, respectively. The B — L number in a comoving 
volume is given as 

riB-L 3 g^ Tr 

= 3 Ci. 3.16 

s 8 g^s f^i 

After the sphaleron process, the B number in a comoving volume is given as 

B = — ~ 0.35 ~ - — ei, (3.17) 

s s orrii 

where we took g^, ~ g^:S- Therefore, the Baryon number is given by 



B = fsmd = ; 7771/ — ^fsmd)~4xlO \ — —V, 3.18 

1287r m, 25QVs:32n^-^gl/^y rn, ' V m/ ' ^ ^ 

where the last similarity is given by taking /sin 5 ~ 0.1 and (7* = (9(100). Roughly, if we 
take 

y^^Mpi/rUi ~ 3 x 10-^ (3.19) 

then we can obtain B ~ lO^^^. The additional condition Fj > {avi)ni becomes 



'^i \ . „ 12 



5 
> 10-1'/ (3.20) 



Mpi 

by using eq. (I3.15p . From the eq. (I3.19p . the additional condition (I3.20p is rewritten as 

mi > 2 X lO-^Mpi ~ 2 X 10^^ GeV. (3.21) 

When rrii ~ 10^^(101^) GeV, the eq. flTTQ]) results in y ~ 2(4) x 10"^ Then F, ~ 2 x 
10-^(7 X 10-^)mi ~ 2 X 10^(7 x 10^) GeV, and therefore, T/j ~ 3 x 10^^(7 x lO^^^ GeV. The 
higher dimensional coupling is given as A/A ~ Ay/Gny/rrii ~ 6 x 10~^^ (10~^'^)(GeV)~^ . 

In this section, we have shown that the Baryon asymmetry in the universe can be 
explained by the B — L production by the decay of some superheavy particle which can 
exist in some GUT models. . 



4 Discussion and Summary 

The initial density of the superheavy fields may be produced non-thermally like the 
preheating [11] and dominate the density of the universe. But here we consider the ther- 
mal abundance of the superheavy fields. If we take (avi) ~ 0.01 — 5-, g^ = g^s = 



mi(GeV) 10^^ 10^^ ^g^^ 10^^ 10^^ lO^*^ Hot relics 


Yi{xlO-^) 0.019 0.10 0.46 1.4 2.0 2.1 2.2 



Table 3: Thermal abundance Yi with y = 0. 



106.75, and g^ff = 1, then the numerical calculation shows that when the mass rrii is 
larger than 10^^ GeV the number density of the i particle behaves like hot relics as Yi ~ 
2 X 10~^ as in Table [31 In the numerical calculation, we used Boltzmann equations with 
Maxwell-Bolzmann approximation for the distribution function. Therefore, for this mass 
ranere, the calculation in case A is reasonable. But if ^^^^^ = |££Zi£±£^y. > i then the 
calculation in case B is preferable because of the entropy production due to the decay of 
the particle i. Therefore, we conclude that thermal abundance of the superheavy fields is 
sufficient to explain the Baryon assymmetry in the universe. One of the point is that since 
the particles are superheavy, the Hubble expansion rate becomes so high that even gauge 
interactions can be out of equilibrium and not affect the generation of asymmetry. This is 
quite different from the usual leptogenesis. 

If one of the right-handed neutrino masses is smaller than the decay temperature and 
10^^ GeV at which the sphaleron process becomes thermalized, then the produced B — 
L number is washed out by the equilibrium of the B — L violating neutrino process and 
the shaleron process[12]. Therefore, the mass of the right-handed neutrinos must be larger 
than the decay temperature in order to obtain the Baryon assymmetry by this mechanism 
if the right-handed neutrino is lighter than 10^^ GeV. 

The B — L violating dimension 7 interactions via the superheavy fields, whose couplings 
are yX/{mfA), can induce the instability of the nucleon. However, the contribution is neg- 
ligible because the effective dimension 6 couplings become very small as y\{hj:,) / {m\h) ^ 

In this paper, we do not introduce the supersymmetry(SUSY), but the extension to the 
SUSY models is straightforward. In some SUSY GUT modelsp3], there may be superheavy 
fields, 10 + 10, 24, and 5 + 5 of SU{5), some of which may produce the Baryon assymmetry 
in the universe, though the serious gravitino problem must be taken into account [H]. 

In this paper, we have studied the possibility that the decay of the superheavy particles, 
which may be induced in grand unified theories as extra fields, produces the non- vanishing 
B — L number, which converts to the Baryon assymmetry by the shaleron process. We 
have shown that if the mass of the superheavy field is larger than 10^'^"^^ GeV, the Baryon 
assymmetry in the universe can be explained by the decay of the superheavy field with 
appropriate couplings. 



Acknowledgments 



We thank M. Tanabashi and T. Yamashita for valuable comments. N.M. is supported 
in part by Grants-in-Aid for Scientific Research from MEXT of Japan. This work was 
partially supported by the Grand-in-Aid for Nagoya University Global COE Program, 
"Quest for Fundamental Principles in the Universe: from Particles to the Solar System 
and the Cosmos" , from the MEXT of Japan. 



A The Calculation of the Mean Net B — L Number 

In this appendix, we will calculate the mean net B — L number e, defined by 

Ci = ^Xi_,/(ri^/ - r-i^j). (A.l) 

/ 



Firstly let us simplify flA.ip . The sum of the branching ratios of a group of decay modes 
g which have the same B — L charges Xi^g = const can be written 

9 fj^9 

Therefore, (lA.ip is transformed to 

ei = Xi^g^{ri^g-rj_,g) + ^x^^f{r,^f-rj^f) (A.3) 

9 fj^g 

= -Xi^g^{ri^f-rj^f) + ^x^^f{ri^f-rj^f) (A.4) 

f7^9 f^9 

= '^ixi^f-Xi'^9)iri^f-ri^f)- (A-5) 

Therefore, in calculating e^, we do not have to calculate all the branching ratios. 

In our case, modes / run two and three body decays which are induced by dim. 4 and 
5 interactions respectively. Here, we choose the three body decays as modes g^ Since 
Xi^2bd. — Xi^3bd. = 2 for all species i as is shown in Table [21 we obtain 

e, = 2 5^ (r.^f - r-^j) (A.6) 

f=2bd. 

= 2 J2 (r^-/ - r^^/)/r^ (a.7) 

f=2bd. 

from f lA.Sp . where Ti^f is the partial decay width and Fj is the total decay width defined 

by 

/ / 

Next, let us calculate difference of partial decay width Fj^j — F^_^j. In case of two 
body decays, the width is given by 

r^^/ = TTr^|M^/l% (A.9) 

IGirrrii 

where Ai is the amplitude. Here we assume that the decay products are massless. Using 
the following Feynman rules 

a (SM) a (SM) 

= -iyiab , i --4^ MSM) - ' A 

h fSM) ^^ hn or /it, 



^Of course, it is the same results if you choose the two body decays as modes g. 



the amplitude can be calculate as follows; 



\M,^ 



'ab\ 



+ 




(A.IO) 



2(Pa ■ Pb) 



Viab 






Kcd ^jdb 
' A A 



J'{mi,mj) + 



m, 



Wiabl 



+ 2ReJ2ylby,c 



j,c,d 



\cd ^jdb 

'■ A A 



J^{mi,mj) + 



(A.ll) 
(A.12) 



where J-" is loop function. On the other hands, the amplitude for the anti-particle can be 
obtained by taking hermite conjugated couplings from the amplitude for the particle; 



\M 



i^ab\ 



'^{Pa-Pb) 



j,c,d 



t , Y^ t ^icd ^jdb -T-i \ , 

Lb + 2^ y]ac^^^^"^i^ ^j) + 



m- 



I |2 I oo \^^ t ■^icd^jdb -r-T/ \ , 

\yiab\ + 2Re 2_^ yUyjac— —-r^ (mi, m^) + 



i,c,d 



A A 



"11 ""J ) 



(A.13) 
(A.14) 



Using flA.Qp . flA.12p and flA.14p . the difference of partial decay width can be written as 

1 



-L 7 — ^nh -L A 



>ab 



IGnrrii 



\M 



i^tab 



\M,^-abV) 



5 E i^ ylbyjc 



Kcd ^jdb 



lm.J^{mi,mj) + 



J,c,( 



Here, ImJ-" is given by 



ImJ-" 



1717 



:f {m]lrr?^i 



(A.15) 
(A.16) 

2(8vr)3^ y-^'-^n (^-17) 

where / is the function defined by fl3.6p . Note that the function T is diverging but IniJ-" 
becomes finite. This is because the imaginary part can be estimated just by tree diagrams 
if Cutkosky rules are applied. 

Finally, we can obtain the mean net B — L number using (lA.Sp . (IA.16P and (lA.lTP : 

'* " WF^^^i ^ ^'^{ylby3ac\cd\]di) f {m]/m]) . (A.18) 

j,a,---,d 
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